Constructing networks of defects with scalar fields 
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We propose a new way to build networks of defects. The idea takes advantage of the deformation 
procedure recently employed to describe defect structures, which we use to construct networks, 
spread from small rudimentary networks that appear in simple models of scalar fields. 
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Networks are of great interest in physics in general. In 
high energy physics, networks appear in diverse contexts, 
usually in scenarios which require the presence of topo- 
logical defects, as junctions of domain walls [l|, cosmic 
strings and branc tiling 0]. The presence and evo- 
lution of domain walls and domain wall networks have 
been investigated in several ways in 0, H, [H, and the 
dynamical evolution of domain wall networks in an ex- 
panding universe has been recently studied in computer 
simulation in Ref. 0]. 

In the present Letter we focus attention on kink net- 
works, that is, we deal with models described by scalar 
fields, which develop spontaneous symmetry breaking of 
discrete symmetry [l|, Q ■ We then take advantage of the 
deformation procedure introduced in Q , and extended to 
other scenarios in Q , to deform a given model, described 
by a potential containing a rudimentary set of minima, 
to get to another model, with the potential giving rise 
to a different set of minima, which may replicate period- 
ically As a bonus, the deformation procedure also gives 
the defect structures of the deformed model in terms of 
the defect solutions of the original model. Thus, in the 
lattice of minima of the deformed model we can nest a 
network of defects in a very natural way. 

This is the main idea underlying this paper, in which 
we use the deformation method to investigate two im- 
portant possibilities, one described by a single real scalar 
field, giving rise to linear networks, and the other by 
a complex scalar field, giving rise to planar networks. 
We focus mainly on the generation of kink-like networks 
described by the deformed models, which are generated 
from simple models, which engender rudimentary net- 
works. 

The idea of constructing networks of defects is not new, 
but the novelty here relies on the use of the deformation 
procedure as a simple and natural way to generate net- 
works. The mechanism is powerful and suggestive, and 
fully motivates the present work. To make the reasoning 
mathematically consistent we consider a model described 
by the Lagrange density with a single real scalar field \ 
in the form 



(1) 



The potential V(x) = (l/2)W' 2 (x) is given in terms of 
the superpotential W — W(x), with the prime stand- 



ing for the derivative with respect to the argument, e.g. 
W'(x) = dW/dx- In this case, the equation of mo- 
tion for static field x — x( x ) carL be reduced to the 
first-order differential equation dx/dx = W'(x)- For 
W = ±(x — X 3 /3) we get the x 4 model, which has the 
set of minima { — 1, 1}. In this case, the defect structure 
represents kink (tanh(x)) or anti-kink (— tanh(x)), with 
energy minimized to the value E = 4/3. For simplicity, 
we are working with dimcnsionlcss fields, space-time co- 
ordinates, mass and coupling constants, with mass and 
coupling constants set to unit. In the one dimensional 
field space, the orbit is a straight line segment which 
connects the two minima. Since the kink or anti-kink 
spans the orbit in the positive or negative sense, they 
may orient the orbit, leading to orientable networks. 

We now use an extension of the deformation procedure 
considered in the first work in [{|. The deformed model 
is described by 

c D = l -d^d»<t ) - l -w' 2 {<p) (2) 

with the deformed potential U(<t>) given in terms of the 
new superpotential W{<f>) = W (f '((j))) / '/' '(</>) . Here /(</>) is 
the deformation function, and we consider f((f>) = sin(0), 
with inverse / T 7 1 (0) = (— l)"Arcsin(</>) 4- nir, with n = 
0, ±1,±2, ... This gives another model, the sine-Gordon 
model with W'(</>) = cos(0). The set of minima is now 
given by {(2n — l)7r/2, (2n + 1)tt/2}. It forms a lattice 
in the entire field space, and n = identifies the cen- 
tral sector with minima {— 7r/2, 7r/2}, n = 1 the sector 
{7r/2, 37r/2}, and n — — 1 the sector {— 37r/2, — 7t/2}, etc. 
The orbit of the original model is now mapped into dis- 
tinct orbits of the new model, giving rise to a specific 
network, which appears as a spreading of the original set 
of two points into the entire field space, the real line in 
the present case. This is illustrated in Fig. 1. 

The above study allows the construction of a regular 
lattice, in which pairs of adjacent minima arc equally 
spaced and connected by kinks and anti-kinks with the 
very same profile and the same energy Ed = 2. We can 
change regularity of the lattice changing the deformation 
function. We take for instance f a (4>) = cos(0 a ), with a 
real and positive, a = 1 leading us to a model similar to 
the former model. It introduces the potential 
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In this case, the set of minima is given by (f> n — ±(nir) 1 / a , 
n = 0, 1, 2, and the distance between consecutive min- 
ima in the lattice increases for a < 1 and decreases for 
a > 1, as we get away from the central minimum at 
the origin. This case gives another tiling, for which the 
distance between minima and the corresponding defect 
energy vary in a nice way, controlled by the parameter a. 
The energy for a = 2/3 in the sector labeled by n is now 
E% 3 = (9/4)(2n+ 1)tt, which increases linearly with n. 
See [l(| for further details. 




Figure 1: Plot of the minima of the original and deformed po- 
tentials, showing how the deformation function projects topo- 
logical sectors in the two models. 

Let us now move to the plane, considering another 
model, described by a single complex held, = 
Xi(x, t) + %X2{x, t), written in terms of the two real holds 
Xi(x,t) and x^ix^t). The specific model which we con- 
sider is described by the Lagrange density 
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(4) 



where the bar stands for complex conjugation. We spec- 
ify the model choosing W^x) as the holomorphic function 



W( X ) = X 
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(5) 



This is the Wess-Zumino model. It was investigated be- 
fore in [ll|, [l2|, EH . The case with N = 3 is interesting 
and illustrative: the vacua manifold has the three points 
Xk = exp(27ri(fc — l)/3), with k = 1,2,3, which depict 
an equilateral triangle in the field plane. And the static 
solutions satisfy the first-order ordinary differential equa- 
tion 



dx 
dx 



-± = e ta W>(x) = e la (l-X%x)) 



(G) 



together with the accompanying complex conjugate, 
where e ta £ S 1 is a phase. We can write W a (x) = 
e~ la W{x) to get d(W a - W a ) = 0. This implies that 
the kink orbits arise when the imaginary part of the su- 
pcrpotcntial is constant 



Im i 



( X 0r)~X 4 (x)/4)) 



const 



(7) 



As the kink orbits connect minima of the potential, 
this constant must also be equal to the value of Im W a 
at those minima, which are the roots of unity. This 
means that sin(2(/s — l)7r/3 — a) = constant. Of course, 
this constant value should be the same at the two dif- 
ferent minima connected by the orbit. Then we have 
lmW aik3) (x {k) ) = lmW a{kj) (x ii} ), and so 



a 



W = -arcsin(cos(7r(fc+j-2)/3)), k > j (8) 



Note that = a (fcj) + tt if j < k. 

We can also use the first-order equations to obtain 



J = \W'(x(x))\ 2 



|(i-x 3 (*))| 



(9) 



where s stands for the "length" on the kink orbits (J7J) - 
see Ref. The kink profiles are then obtained by in- 
verting these relations between the real part of the super- 
potential and s. The energy of the static configurations 
is £7 =(3/2) |sin((fc-j>/3)|. 

We now turn attention to the deformation procedure. 
We follow the second work in 0. It is interesting to 
express the deformed system in terms of another complex 
field, 4>{x,t) = <fii(x,t) + i<j)2(x,t), related to the original 
one by means of a holomorphic function / = f(4>) such 
that 

X = f{4>) = fi(4>i,h)+if2^t,h) (10a) 
dh = dh, dh^ dh (1Qb) 
d<t>i dfy-i d<j)2 d4>i 

The deformed Lagrange density has the form 



c D =h^-^m (id 

2 ' /W(0 

The deformed model Ld can be defined by the new su- 
pcrpotcntial 



W{4>) 
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[12) 



In this case, the "deformed" first-order equations are 



-^=e M W'(0); / = 
ax dx 



e- ia W'{<t)) (13) 



The defect solutions for this system are obtained from 
the solutions of © by simply taking the inverse of the 
deformation function: <p K (x) = f^ 1 (x K { x ))- Thus, we 
can establish the following relation between the deformed 
and original equations: if X ( x ) 1S a kink- like solution of 
the original model, we have that 

lm.W{x K {x)) = const; ReW( X K (x)) = s (14) 

and then 4> K (x) = f' 1 (x K ( x )) is kink- like solution of the 
deformed model, obeying 



ImW(/- 1 ( X K (x))) = const; ReW(/" 1 (x K (»))) =<r 



where a is defined by 



*= \ \w{r 1 {x K {x)))\ 2 dx 



(15) 



(16) 



Although the method is general, we now specify the 
deformed model choosing f(<f>) = W(<fr). This constrains 
the function /(</>) to obey the equation 



/ww = (17) 
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A function / satisfying this condition provides a potential 
U((j>, cj>) for the new model which is well defined (finite) 
at the critical points of /(</>), e.g. the zeros of f'(4>). As 
a bonus, the procedure leads to a very simple expression 
for the deformed superpotential. We change x f° r /(</>) 
in the general expression ([5|) to get the potential 
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v = -(i-.r(<t>m-f N m 

As stated in p7|) . we can then choose 

/' 2 (0) = (-l) N (l-/ N (0)) 



(18) 



(19) 



The solution of this equation solves the general problem. 

We illustrate the general results with N = 3. Here we 
have 



f' 2 ^) = / 3 (0) - 1 



(20) 



The solution is the equianharmonic case of the Weier- 
strass V function 



W) =/(</>) =4^(4-30; 0,1) 



(21) 



The Weierstrass V function is defined as the solution of 
the ODE 



{V'{z)f = 4V 3 (z) - g 2 V(z) 



93 



(22) 



The elliptic function which solves the differential equa- 
tion above is doubly periodic function defined as the se- 
ries 



V(z) 
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(z — A(m,n)) 2 A(m,ny 



(23) 



where A(m, n) — 2m,LO\ + 2nt03, with m, n € Z and 
m 2 + n 2 7^ 0. Therefore, the deformation function is, up 
to a factor, the Weierstrass V function with invariants 
g-2 = and g^ — 1, and we denote it by Vqi(z). This 
function is meromorphic, with an infinite number of poles 
congruent to the irreducible pole of order two in the fun- 
damental period parallelogram (FPP). 

Here we get W a {4>) = e _M 4^ Vm(^~ 3 <f>), and so the 
deformed potential can be written as 



tf(&0) = ^oi(4"fy)Poi'(4- 



(24) 



The potential spans the plane replicating the triangular 
structure as shown in Fig. 2, and in Fig. 3. 

The new potential is doubly periodic with an structure 
inherited from the "half-periods" of V . The set of zeros 
of the potential in the FPP has three elements i^ 1 ) = 
wi = wa(l/2 - »V3/2), <t> (2) = w 3 = wa(l/2 + iy/3/2), 
and (j)^ =lo 2 = 4^r 3 (l/3)/47r. The set of all the zeros 
of U form a lattice which tile the entire field plane, as we 
show in Fig. 3. 

The potential obtained from the deformation proce- 
dure has the same number of zeros in the FPP as the 




Figure 2: The case N = 3, showing the potential U{4>,(j>) 
(upper panel) and its mechanical analogue —U(<j>, <j>) near a 
pole (lower panel). Note that in the lower panel the zeros are 
now maxima. 



original model in the whole field space. Besides, one pole 
of sixth order arises at the origin due to the meromor- 
phic structure of 7 3 q 1 (4 _ 3 <fi); see Fig. 3. However, this 
structure is infinitely repeated in the deformed model, 
according to the two periods lo\ and UJ3 determining 
the modular parameter r = UJ3/0J1 = — 1/2 + iy/3/2 of 
the Riemann surface of genus 1 associated with this V- 
Weierstrass function. 

We now compare the 'P-kink orbits with the orbits of 
the original model. If \ K i x ) IS a solution of and © 
then <j) K (x) = 4 1 V^ 1 (4" i% K 0*0) s °lves 



Ime- M 4^oi(4"^ K (x)) 
Re 4^oi(4-^(:e)) 



where 



7> 01 (4-^(*)) 



const 
a 



dx 



(25) 



(26) 



Since the deformation function is a conformal transfor- 
mation, angles arc preserved, and the same values of a 
as in the non deformed case give the kink orbits. There 
are three types, and here we just inform that the nearest 
neighbor type (12), (23) and (31) minima are connected 
by orbits which follow specific sequences - see Ref. [lij ]. 
We notice that the kink orbits go around, circumvent- 
ing the singularities which stand at the center of circles 
depicted by the orbits themselves. Like in the case of a 
single real scalar field, we can also break the lattice reg- 
ularity in this case - see [l4| for further details on this 
issue. 

We further illustrate the problem with N = 4. Here we 
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Figure 3: The case N = 3, showing the minima (•) and poles 
(o) of the potential (upper panel, left) and the kink orbits 
in the FPP (upper panel, right). The lower panel shows the 
lattice of minima and poles of the deformed potential and the 
accompanying network of kink orbits. 

have W(x) = X ~ X 5 /5- Thus, the potential is 

V(X,X) = 1(1-X 4 )(1-X 4 ) (27) 

We write x = /(0) to show that the special deformation 
function must satisfy 

fW/1 = V(W 4 (0)(1-7*W)) (28) 

As before, we choose the holomorphic solution of 

f 2 (0) = l-/ 4 (0) (29) 

The solution is the elliptic sine of parameter k 2 = — 1, 
the Gauss's sinus lemniscaticus /(</>) = sn(c/>, — 1). The 
new superpotential is W a (4>) = e~ t<x sn{4>, — 1) and the 
deformed potential then reads 

U{4>,4>) = ||cn(0,-l)| 2 ■ |dn(0,-l)| 2 (30) 

The new potential is doubly periodic with an structure 
inherited from the "quarter-periods" K (— 1) = Wi/4 and 
iK(2) = W2/4 of the twelve Jacobi elliptic functions. 
Here K(—l) rs 1.31103 is the complete elliptic inte- 
gral of the first type, a quarter of the length of the 
lemniscate curve in field space: (</> 2 + 4>?,) 2 = — </>§. 
K(2) R 1.31103 — il. 31103 is the complementary com- 
plete elliptic integral of K(— 1). 

The set of zeros of the potential in the FPP are 
<f>W = wi/4,0( 2 ) = iujjA,^ = -cji/4,<^ 4 ) = -iwi/4, 
whereas the set of all the zeros of U form a quadrangu- 
lar lattice in the whole configuration space. This is de- 
picted in Fig. 4 and will be fully considered in Ref. [lil ]. 



Differently from the former case, however, here the or- 
bits may connect the minima in two distinct ways: one, 
with curved lines, in the sequence (1, 2), (2, 3), (3, 4), and 
(4, 1), and the other with straight line segments, in the 
sequence (1,3) and (2,4) - see [lj] for further details on 
this issue. 




Figure 4: The case N = 4, showing the set of minima (•) and 
poles (o) of the deformed potential and the accompanying 
network of kink orbits. 

In summary, in this work we have used the procedure 
developed in |§, 0] to deform a given model in a way such 
that its set of minima could be replicated in the entire 
field space. The idea was developed in the real line, for 
the case of a real field, and in the plane, for the case of 
a complex field. Since the set of minima are connected 
by algebraic orbits describing defect structures in field 
space, we have also been able to replicate the algebraic 
orbits in the entire field space of the deformed model, 
naturally building networks of defects, which are spread 
from rudimentary networks into the entire field space. 

A natural extension of this work concerns the construc- 
tion of irregular lattices and networks in the plane, in the 
case of a complex field, which we will study in our next 
work, now under preparation [l4j|. Another extension 
concerns the use of three real fields, to investigate if it is 
possible to tile the space in a way similar to the case of 
planar networks here considered. 

We recall that a kink-like defect in general splits the 
space into two distinct regions, so we could also think as 
in [l[, using two spatial dimensions, to see how the kinks 
orbits that we have just obtained could tile the plane 
with regular and/or irregular polygons, with triple junc- 
tions for N = 3, and with quartic junctions for N = 4. 
Another interesting issue could address the same prob- 
lem, but now embedding the scalar fields in a curved 
space-time, following the lines of Ref. [J]. This would 
lead us to another route, in which we could try to under- 
stand how the networks here introduced would change in 
a curved background. We can also think of making the 
space-time dynamically curved, to see how the domain 
wall networks could follow the evolution investigated in 
0- These and other related issues are presently under 
consideration, and we hope to report on the them in the 
near future. 
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